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Abstract 

The multiquark confining potential is proportional to the total distance of the fun- 
damental strings linking the quarks and antiquarks. We address the computation of 
the total string distance an of the Fermat points where the different strings meet. 
For a meson (quark-antiquark system) the distance is trivially the quark-antiquark 
distance. For a baryon (three quark system) the problem was solved geometrically 
from the onset, b)' Fermat and by Torricelli. The geometrical solution can be deter- 
mined just with a rule and a compass, but translation of the geometrical solution to 
an analytical expression is not as trivial. For tetraquarks, pentaquarks, hexaquarks, 
etc, the geometrical solution is much more complicated. Here we provide an iterative 
method, converging fast to the correct Fermat points and the total distances, rele- 
vant for the multiquark potentials. We also review briefly the geometrical methods 
leading to the Fermat points and to the total distances. 
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1 Introduction 

Fermat proposed to Torricelli the problem of finding the point in a triangle 
minimizing the sum of the distances to the three respective vertices. This first 
Fermat point or Torricelli point [1,2,3,4,5,6,7,8,9], is the isogonic point, since 
in a sufficiently acute triangle the angle formed by the segments connecting 
any two vertices with it is 120 degrees. 

Lately this problem became relevant for quark physics because the multiquark 
confining potential is proportional to the total distance of the fundamental 
strings finking the quarks and antiquarks. Of course, we address here the case 
where we have a single multiquark and not many free or molecular mesons 
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and baryons, where the confining potential would be different. The three-body 
star-like potential has already been used long ago in Baryons [10], however for 
many years there was a debate in the lattice QCD community on the two-body 
versus three body nature of the confining potential for baryons. Recently, the 
study of fiux tubes in Lattice QCD for Baryons (triquarks) by Takahashi et 
al [11] confirmed the three-body star-like confining potential . Very recently, 
the Wilson loop technique was applied to tetraquarks [12] and pentaquarks 
[13] by Okiharu et al, showing that the confining potential is provided by a 
fundamental string linking all the quarks and antiquarks. Cardoso et al also 
confirmed this result with the Wilson loops for hybrids [14] and for three gluon 
glueballs [15]. Thus we assume that the confining component of the multiquark 
potential is. 



where a is the string tension, is the position of the quark or antiquark Qi , 
Ta is the position of the Fermat point F^, and we use respectively arab digits 
i = 1, 2, 3 • • • for the quarks (antiquarks) and roman digits a = I, II, III ■ ■ ■ 
for the Fermat points. Thus the Fermat problem of finding the paths mini- 
mizing the total distance is equivalent to the physics problem of computing 
the multiquark potential. Notice that there are already some proposed exper- 
imental signals of tetraquarks, and the next generation of Hadronic Detectors 
may eventually observe multiquark hadrons. 

The geometries of the strings of the first five multiquarks are depicted in Fig. 
1. Eq. (1) and Fig. 1 extend the definition of the Fermat point of a triangle 
to the Fermat point of polygons in three dimensions with more points. With 
the present definition, where confinement is produced by fundamental strings, 
the strings meet in internal three-string vertices. The number of quarks can 
always be increased replacing a quark (antiquark) by a Fermat point and a 
diquark (di-antiquark). Thus the number of quarks (and antiquarks) minus the 
number of Fermat points is a constant. Since in the meson and baryon this 
constant is 2, the number of Fermat points is — 2 where N is the number of 
quarks and antiquarks. Moreover in eq. (1) we are only suming over distances 
between points linked by strings. 

For a meson (quark-antiquark system) the distance is trivially the quark- 
antiquark distance. For a baryon (three quark system) the problem was first 
solved geometrically by Fermat and by Torricelli. In the case of 3 quarks, the 
minimization of the potential in eq. (1) implies that. 



(1) 



rii + r2i + r3i = , 



(2) 



2 




Figure 1. The geometries of the string sections linking the first five multiquarks. 
Notice that the number of Fermat points Fa is — 2 where N is the number of 
quarks and antiquarks Qi in the multiquark. 

and it is clear that the solution is that, either the triangle is not sufficiently 
acute, or the angles are all equal to 120°, 

rTiTr^i = rTzTrs/ = r^JTrT/ = 120° . (3) 

Due to the beauty of the triangles, and also to their simphcity, there are 
numerous geometry textbooks and articles on the Fermat - Torricelli point 
[1,2,3,4,5,6,7,8,9]. However, when the number of quarks increase, to tetraquarks, 
pentaquarks, etc, the geometric construction of the Fermat points becomes 
more and more difficult. Thus a numerical solution of this problem is wel- 
come. 

Here we address the computation of the total string distance and of the Fermat 
points where the different strings meet. In Section 2 we review briefly the 
geometrical methods leading to the Fermat points and to the total distances. 
In Section 3 we provide an iterative method, converging fast to the correct 
Fermat points and the total distances, relevant for the multiquark potentials. 
We detail the cases of the baryon, the tetraquark, the pentaquark and the 
hexaquark. In Section 4 we conclude. 
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Figure 2. A step in the geometric construction of the first Fermat point of an acute 
triangle. Starting from the segment Q\ Q2, an equilateral triangle with vertex V12 is 
constructed. The Format point Fj belongs to the arc of circle centered in V12 and 
passing by Qi and Q2- 

2 Brief review of the geometrical method 

In an acute triangle, the Fermat point Fj is the isogenic point, defined in eq. 
(3). To construct the isogonic point, we start by the first pair of vertices Qi 
and (52, noticing that the set of points Fj with fixed angle Q1F1Q2 = 120° 
belong to an arc of circle. Moreover this circle is centred in the other vertex V12 
of an equilateral triangle including Qi and Q2 ■ In Fig. 2 we show the 120° arc 
of circle, the equilateral triangle, and a segment including the points V12 and 
Fj. Notice that the other end of this segment forms with the segments Qi Fj 
and Fj Q2 angles of 120°. Thus the isogonic point belongs this arc of circle. 
This point is at the intersection of the segments V12Q3, ^sOi and V?nQ2- The 
construction of first Fermat point Fj is illustrated in Fig. 3. It is very simple 
in a geometrical perspective, and it can be done just with a compass and a 
rule. 

We now proceed with the tctraquark. This geometrical method can be ex- 
tended to construct the two Fermat points Fj and Fjj of a tetraquark. Notice 
that in the tetraquark we have four points, and thus in general the points 
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Figure 3. The geometrical method to construct the Fermat Fj point of a Baryon. 
The Fermat point is at the intersection fo the three segments Vi2Q^, V2zQi and 

Qi, Q2, Q3 and (^4 are not coplanar. Thus the vertices V12 and V34 are not, 
from the onset determined, only the circles where they belong are determined 
with the technique already used for the baryon. To determine the vertices, 
notice that the vertex V12 must be as far as possible from the segment Q3 Q4 
and that the vertex V34 must be as far as possible from the segment QiQ2- 
Thus we find that the segment V12 V34 must intersect the segment Qi Q2 and 
the segment Q3Q4. Then, once the segment V12 V34 is determined, the Fermat 
points Fj and Fn are determined because the distances V12 Fi — Qi Q2 and 
V34F// = Q3Q4. This is illustrated in Fig. 4. 

Although the solutions are simple geometrically, the algebraic computation 
of the cartesian coordinates of the Fermat points Fa is not simple. Let us 
consider the baryon case of triangle in a three dimensional space. We first 
need to check whether the triangle is acute. If the triangle is acute, we find 
each of the vertices Vij of the equilateral triangles with a system of three 
equations, two hnear equations stating that the vertex belongs to the plane 
of Qi, Q2 and Q3 and to the plane of the mediatrices of Qi and Qj , and one 
quadratic equation to fix the distance of the vertex to the medium point of 
Qi and Qj. We have to ensure that the vertices Vij point outwards the initial 
triangle (5i<52<53- To find the Fermat point Fj we need to get two vertices. 
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Figure 4. The geometrical method to construct the two Fermat points Fi and Fn of 
a tetraquark. Notice that The points Qi, Q2, Q3 and Q4 are not coplanar. Thus the 
vertices V12 and Vm are not, from the onset determined, we only the circles where 
they belong. Nevertheless it is possible to determine tehm geometrically, knowing 
that the segment V12 V34 intersects the segment Qi Q2 and the segment Q3 Q4. Then, 
once the segment V12 V34 is determined, the Fermat points Fj and Fn are determined 
because the distances V12 Fj = Q\ Q2 and V34 Fu = Q3 Qi- 

say V12 and 1^23, and then to intersect the segments Vi2Q3 and V23Q1. So we 
first need if statements, and in the acute case we have to solve a total of two 
quadratic equations and seven linear ones. 

While this system of equations and inequations is exactly solvable for a tri- 
angle, the algebraic method gets quite difficult for larger multiquarks. Thus 
another method would be welcome to compute the Fermat points and the 
total distance d — J2i,a''^ia- 



3 Iterative method 



We now propose an iterative method, designed for a computational determi- 
nation of the confining potential of multiquarks. 
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Figure 5. Convergence of the numerical iterative method method to construct the 
two Fermat points F/ and Fjj of an arbitrary tetraquark. The quarks are depicted as 
circles, the antiquarks as squares, and the converging Fermat Points as diamonds. 
We show projections in the xy, yz and zx planes. Visually, after 6 iterations the 
iteration has converged. On average, a 10~^ precision is achieved for the tetraquark 
after 20 iterations. 

We start from the triquark case of the baryon. We use the notation, 



Minimizing the total distance ru + r2i + r^j with regards to the three coor- 
dinates of the Fermat point F/, we get the system equations, 




{xi,yi,zi) , 



{xi - XiY + {yi - yif + {zi - ZiY . 
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which is non-hnear. To solve these algebraic equations, of the form 



X = X{x,y,z) , 
y = Y{x,y,z) , 

z = Z{x,y,z) , (6) 



an iterative method can be used. With a relaxation coefficient cu, we iterate 
the series, 



Xn+l^UjX{Xn,yn,Zn) + {I - U;)Xn , 
yn+l^U;Y{Xn, yn: Zn) + (1 " ^)yn , 

Zn+l^UjZ{Xn, yn, ^n) + (1 " ^)^n , (7) 

starting with, as an initial guess, the barycentre 



ri + r2 + ra 

r/o = o ■ (8) 



A first numerical check shows that the method converges rapidly to the Fermat 
path of the triangle. We get get results accurate to a precision of 10~^, for the 
total distance d — ru + r2i + r^i after a number of iterations of the order of 
12, for an optimized relaxation factor cu — 1.7. 

Thus we extend our iterarive method to the cases of the next multiquarks. We 
simply get a system of there equations per Fermat point. For the tetraquark 
the equations are. 
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and we also check that the iterations converge fast to the Fermat paths for the 
tetraquark. We illustrate graphically, for an arbitrary geometric configuration, 
the convergence in the case of a tetraquark in Fig. 5. For the pentaquark the 
equations are. 
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for the hexaquark the equations are. 
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Tabic 1 

Convergence of the iterative method for different multiquarks, based on a million of 
random generated quark positions, for each multiquark, and for different precisions 
p = where d is the total Fermat distance. The values of the relaxation factor a; 
are the ones which minimize the convergence time up to a precision p = 10"^. 
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Figure 6. Distribution of the number of iterations needed to converge to a preci- 
sion p = 10~^, based on a million of random generated quark positions, for each 
multiquark. 

and for larger multiquarks the extension of these equations is obvious. 

The convergence of the method for the first multiquark systems is shown in 
Table 1 and in the Fig. 6 

We study numerically the convergence of the method, with a randomly gener- 
ated sample of 10^ geometric configurations, for each of the multiquarks. We 
first optimize the relaxation factor uj in order to reduce the needed number 
of iteration steps, to get results accurate to a precision of 10~^, for the total 
distance d = Y.i^a'^ia- Since the number of iterations is not a constant, we 
demand the uj minimizing the average d over the sample of geometric config- 
urations. The optimized uo and the resulting d are shown in Table 1 for the 
baryon, the tetraquark, the pentaquark and the hexaquark. Notice that the 
convergence is fast, even for the hexaquark. 

Then we compute the distribution of the number of convergences we have per 
number of iterations. The number of configurations is large, and so we don't 
need to join different numbers of iterations in a bin. In Fig. 6 we show the 
distribution of the number of distribution of the number of convergences we 
have per number of iterations for the baryon, the tetraquark, the pentaquark 
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and the hexaquark. 



4 Conclusion 

We study an iterative method to find the Fermat points and Fermat distances 
in multiquarks. This method replaces the geometrical method, which has only 
been applied so far to the baryon (triquark). 

The method is very simple to programme, and it converges both in the case of 
acute angles (smaller than 120°) and of larger angles. This avoids the problem 
of checking for non-acute angle, simple for a Baryon, but harder for the larger 
multiquarks. 

This method is suited to be used in the solution of the Schrodinger equation 
in a quark model, where we have to compute the potential for several different 
positions of the quarks and antiquarks. 

Thus we use, as a convergence criterion, the precision of 10~^, one part per 
million, in the total distance of the Fermat paths defined in Fig. 1. Even for this 
extremely fine precision, far beyond the normal precision of quark models, the 
method is quite fast, as show in Table 1. The number of necessary iterations n 
to achieve a precision p = grows quadratically with the number of desired 
correct decimal cases, — logp is proportional to n^. 

The present computational technique enables the precision studies of the mul- 
tiquarks with string confinement in the quark model. 
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